An Estimate of the Branching Fraction of r 
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We calculate the expected branching fraction of the second-class-current decay r — + ng'u T , moti- 
vated by a a recent experimental upper-limit determination of this quantity. The largest contribution 
to the branching fraction is due to the intermediate ao(980) scalar meson, assuming it is a ud state. 
Smaller contributions arise from ao(1450), p(770), and p(1450). Our calculated values are substan- 
tially below the experimental upper limit, and are smaller still if the ao(980) is a four-quark state, 
as often suggested. Thus, a precise measurement or tight upper limit has the potential to determine 
the nature of the ao(980), as well as search for new scalar interactions. 

PACS numbers: 



I. INTRODUCTION 

In a recent paper [l[ , we considered the branching frac- 
tion of the isospin- violating decay t — > -nr]v r . We found 
an expected branching fraction of 



B = B{t -> 7rr]U T ) = (0.3 - 1.0) x 10" 



(1) 



in rough agreement with a detailed chiral-perturbation- 
theory calculation || and other evaluations which 
yielded central values in the range 



B= (1.2-1.6) x 10" 



(2) 



The experimental bound on this branching fraction, B < 
1.4 x 10~ 4 4], was obtained by CLEO with an e + e~- 
collision data sample of 3.5 fb -1 , a fraction of a percent 
of currently available integrated luminosity. The only re- 
lated high-luminosity measurement is a stringent BABAR 
upper limit on the branching fraction of r — > ixrj'y T ||, 



B' = B(t -» mfv T ) < 7.2 x 10~ 6 @90% CL, 



(3) 



obtained with an integrated luminosity of 384 fb^ 1 . 

The fact that the experimental limit is lower than the 
results summarized in Eq. @ raises the question of a pos- 
sible discrepancy between theory and experiment. There- 
fore, our goal in this article is to calculate the expected 
value of B' and compare it to the experimental limit. 
We adapt the methods used in Ref. 1] to the present 
case, noting that a chiral-perturbation-theory calculation 
of this process, as performed for t — > Trr\v T by Neufeld 
and Rupertsberger [2J, would be very useful. 

First, we note several similarities and differences be- 
tween the calculations of B' and B: 

• The uu + dd fraction of the wave function which, 
unlike the ss and gg parts, contributes to the de- 
cay amplitude, may be smaller for the rj'. While it 
appears that the magnitude of the ss part in rela- 
tion to that of the light quarks is very similar for 
both states, the current estimate of the gg fraction 



of the wave function, Z nn . is 



\z gg ? 



0.3 ±0.2 



In our calculations we take Z gg = 0, as this yields 
the most conservative limits on B' ', and since the 
modification for finite values of Z gg is straightfor- 
ward. 

• Calculations of B in Refs. [3, [3,3 rely on extrapola- 
tions utilizing intermediate, low-mass 

and ++ hadrons. Obvious intermediate states for 
the decay r — > ixr\v T are the ground-state mesons 
p(770) and ao(980). In the case of r — > nrj'v T , 
these are off-shell processes, and the contributions 
of these resonances are suppressed. On the other 
hand, we do have now on-shell decays involving 
the next 1~~ and ++ states. These are the 
p' = p(1450) and ao(1450), which contribute to the 
P- and S*-wave components of the decay, respec- 
tively. 

• The p and p' vectors are the quark-model ud, 
S- wave 1 ground state and first radial excita- 
tion, respectively. However, the theoretical assign- 
ment of the ao(980) (and, consequently, that of 
the ao(1450) as well) is ambiguous, generating the 
largest uncertainty in both B and B 1 . Conversely, 
information on these branching fractions can help 
resolve the longstanding dilemma of the "KK- 
threshold" state do (980). The significant branch- 
ing fractions of do (980) and /o(980) decays to KK, 
despite the very small phase space, seem inconsis- 
tent with these mesons being the ground states of 
the quark- model scalar nonet, motivating a four- 
quark (udss) interpretation [7j. In this case, the 
ud scalar ground state should most likely be identi- 
fied with ao(1450). However, this would make the 
scalar 190 MeV heavier than the axial vector state 
cii(1260), implying a pattern oiL-S splitting differ- 
ent from what is observed in any other L = 1, qq' 
system. The more appealing possibility, namely, 
that the two 980-MeV states arc indeed just ud 
states, may have been partially resurrected in re- 
cent work [8|], in which 'tHooft's uuddss six-quark 
vertex was utilized to admix the 2- and 4- quark 
states. 
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The plan of this note is as follows. As we did in 
Ref. [l|, we discuss separately our estimates of the P- 
and S- wave contributions to B' . In Sec HI] we present 
the more robust results for the P-wave part, calculating 
upper bounds on the contributions of the p and p' using 
recently published experimental data involving 77' and r 
decays. In Sec IIIII we present the less clear-cut estimate 
of the S'-wave component. This contribution depends 
most strongly on whether the ao(980) is a 4-quark state 
or the ud ground state. In any event, our predictions for 
B(r — > 7T7/V T ) lie significantly below the BABAR limit ||. 
A brief summary and future outlook are given in Sec lIVI 



II. THE L = 1 CONTRIBUTION 

In Ref. [H, we obtained the L = 1 contribution to B 
assuming that it was dominated by the p, an assumption 
justified by the large branching fraction B(t — > pv T ). We 
compared this branching fraction to B using the ratio of 
coupling constants g np ^ / 'g p7r ir, where <7 pvrTr was related to 
the width of the p, and g W7r was obtained by analyzing 
the Dalitz-plot distribution of the decay 77 — > ir + ir~ir°, 
taking the scalar contribution to 77 — ► 7r + 7r~7r from 
B(j]^ 7r 7r°7r°). 

This procedure is not directly applicable to B' , since 
there is no experimental information on the Dalitz-plot 
distribution of the decay 77' — > ir + ir~ir°, nor a measure- 
ment of #(77' — * 7r°7r°7r ). Therefore, we make use of 
the fact that the branching fraction B(rj' — > 7r + 7r~7r°) de- 



pends on the coupling constant g^pm under the conser- 
vative assumption that the states dominate the de- 
cay 77' — > 7T + 7r _ 7r°. This will yield a conservative upper 
bound on g^pm from which we obtain an upper bound 
on the p contribution to r — > nrfis T . We discuss the 
likelihood of this assumption and its implications below. 

The differential branching fraction of 77' — * 7r + 7r~7r° as 
a function of the Dalitz-plot position is given by 



dV 7 y — > 7r + 7r (^77' p-ir9 p7T7r) Q 



384\/37i 



m 



,r , 



■\M\ 2 dXdY, (4) 



where 



3777-7, 



(5) 



is the kinetic energy in the decay, and 



X 



V3 

Q~ (T+ 



Y ^ Q T ° ~ 1 



(G) 



are the Dalitz-plot variables, with T c being the kinetic 
energy of the pion with charge c. Assuming p dominance, 
we obtain from Eq. (15) of Ref. the reduced matrix 
element 



M 



rY 



\Y 2 + X 2 ) 



ir2(iy2_ x2) - 



(7) 



where 



m V 'Q 



3"V 



iY p mp 



= 1.6 + 0.7?:. 



(8) 



The product (g v 'pirgpTrir) 2 is then found by integrating 
Eq. ([4]) over the Dalitz plot. In the 77 — > 7r + 7r~7r° case, 
we exploited the small value of r to simplify the expres- 
sion by expanding in r. Due to the O(l) value of r for 
7/' — > 7T + 7r~7r , we resort to numerical integration, which 
yields 



|M| 2 dAdy = 2.4. 



(9) 



From this we obtain, using B{rj - 
lO" 3 9] and 9fm7r =6.0 |, 

5^ < 0.025. 



TT+TT-TT ) = 3.7 X 



(10) 



As a cross check, we apply the procedure to the de- 
cay 77 — * 7T + 7r~7r°, obtaining g^p-^ < 0.52. This value 
is to be compared to the one obtained from the more 
precise Dalitz-plot analysis in Ref. [![, g-qpn w 0.085. 
The factor of 6 ratio between the results reflects the 
fact that the procedure used here yields but a conser- 
vative upper bound, obtained by assuming that the de- 
cay 77' — > 7T + 7r~7r is dominated by the p resonances. 
This assumption is manifestly false, as the 77' — ► 7r + 7r~7r 
Dalitz-plot distribution is in much better agreement with 
a flat distribution than with that expected from p^ dom- 
inance Q. By contrast, in Ref. [l|, the value of gf^ p7r 
obtained from the Dalitz-plot distribution yielded good 
agreement between the expected and measured values of 

£(77^ TT+TT-TT ). 

With this point in mind, we proceed to use the up- 
per bound on g^'pn to calculate the upper bound on the 
p contribution to B(t — ► nrfv T ). We do this by relating 
B(t — ► p(- v ,\v T ) to B(t — + Plnri)^) via the ratio of cou- 
pling constants and phase-space factors 



(11) 



where pt V7) n indicates that the p is observed in the 7T?7' 
final state, and V(X) is the integral over the Dalitz plot 
of the three-body decay X. The ratio of phase-space 
integrals is 0.06, with up to 15% variation depending 
on whether one uses Blatt-Weisskopf and s— dependent 
widths for the p and on the choice of angular distribution. 
Using B(t — > P(-Kn) v T) — 3.6 x 10~ 6 [l|, we obtain 



#(r -> p<srn>)V T ) < 2 x 10" 



(12) 



more than two orders of magnitude below the BABAR up- 
per limit, Eq. ©. 

Next, we evaluate the contribution of the on-shell p' . 
One expects that this state, being a radial excitation 
and hence having a node in its wave-function, couples 
to the ground-state particles 77 and tt more weakly than 
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the p. We hypothesize that this p' suppression mecha- 
nism works equally strongly for the final states irq' and 
7T7T, leading to an equality of the ratios of the squared 
matrix elements 

B(t -> p {mi >)V T ) V(t p\ nr , t) v T ) 
B{r -> p\^)V T ) V(t -> p { ^)V T ) 



B{t^p^)V t ) V(t -> p' {7rn) u T )' 
The relevant phase-space integral ratios are 



(13) 



^0" P\ ml ') v T) 

V(t ^ P{^)Vr) 
V(T - p' ( ^ T ) 



0.06, 



2.5. 



(14) 



We use the upper bound of Eq. (TT2"]) and the central value 
plus one standard deviation of the recent Belle result (lpj 



= 0.15 ±0.05 



y B(r -> P{iw)V r ) 
to obtain the conservative upper limit 

B(t ^ p 1 ( ^, )Vt ) <8x 10- 8 



+0.15 
-0.04 



(15) 



(16) 



We note that this is an upper bound both due to the way 
we use Eq. (fT5|) and since Eq. (TT2"]) is an upper bound. 



III. THE L = CONTRIBUTION 

Calculating the L = contributions to B' is not as 
straightforward as the L = 1 case, where one can make 
use of the dominant p coupling to the leptonic vector 
current. Therefore, is is important to evaluate the scalar 
component using different methods, as has been done for 
the r — > ~kt]v t decay 0,0,13. It should be noted that 
these calculation are performed under the assumption 
that the relevant scalar resonances are ud states. The 
coupling of a 4-quark state to the ud scalar current is 
"Zweig-Rule" suppressed, making it significantly smaller 
than the predictions. 

Here we perform a more detailed version of the calcula- 
tion used in Ref. [l| . We begin with the ratio of branching 
fractions 



R a ° 



B(t — > a v T ) 
B(t — ► a\v T ) 

Pa „ \(aa\V h r\0) (v r \J? 



P«i \{ ai \A hlJ \Q) {v t \J?\t)\' 



(17) 



where ao stands for either ao(980) or ao(1450), ai is the 
ai(1260), px is the T-rest-frame momentum of the prod- 
ucts of the decay r — ► Xv T , = r ipu(x)lij,'4'd(x) is the 
hadronic vector current, = 4i u (x)l[i'y b '4'd(x) is the 



hadronic axial vector current, and = ■ip UT (x)j fJ '(l — 
j 5 )ip T (x) is the leptonic current. The calculation of the 
leptonic parts of this ratio is well defined, while all the 
uncertainty in the hadronic parts comes down to a single 
parameter £, which shall be defined shortly. With this in 
mind, we can take the ao matrix element to be 



{a \V hfl \Q) = f -^ {a o \S h \0) , 



(18) 



where /o is an isospin- violation suppression factor, and 
Sh = ip u (x)ipd(x) is the scalar current operator. The 
weak vector current is conserved up to the difference be- 
tween the u- and d-quark masses, plus a smaller electro- 
magnetic part that we neglect. Therefore, 



<9 M ^V w (m d ~ m u )S h . 
Using this relation in Eq. (TTSJ) yields 



fo = 



m d 



(19) 



(20) 



We use the fact that both the ao and the a\ are P-wave 
states to relate the axial and scalar decay constants 



<ai|A„|0>=£e* (a |5|0) 



(21) 



We note that this is reminiscent of applying SU(6) [Tl| 
or, in this case, just SU(4) [lj] flavor-spin symmetry to 
the (L = 0) 15-plet plus singlet containing the w, p, 7], 
and ui, or the (L — 1) states ao, Oi, fo, and h\. 

Naively, one expects £ in Eq. (f2TTl to be of order 
unity. However, this parameter incorporates all the 
hadronic uncertainty in our procedure. With Eqs. (1181 
[2Tj) . Eq. (TIT)) becomes, after spin averaging and index 
contraction, 



Pat V m 



m„ 



a a 

m 



ai 



III,, 



a J I + 2 (m T / m ai Y 



(22) 



This yields the branching fractions 



B(t -► a (980)z/ T ) = 1.6 x 10~ 6 |£| 2 , 
B(t -> a (1450K) = 6.4 x 10~ 8 |£| 2 , (23) 

where, as in Ref. [l|, we chose the mass difference of the 
two light quarks to be 4 MeV [l3| and, assuming that 
the t — > 3iri> T decay is dominated by the oi, we took 
B(t -> oii/ r ) = 0.18. We compare B{r -> a (980)i/ T ) of 
Eq. (gU with the value B = 1.2 x 10" 5 obtained from 
the more elaborate calculation of Ref. [31, minus the p 
contribution to B, which is 3.6 x 10~ 6 [l|. This yields 
|£| 2 ~ 5, from which we conclude 



B(t -> a (U50)v T ) w 3 x 10~ 7 . 



(24) 



The ao(1450) contribution to r — > -kt{v t depends also 
on the branching fraction S(ao(1450) — > 7177'), regarding 
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which there is only partial information. However, from 
the branching-fraction measurements that have been 
made it is clear that £(a (1450) -> irr)') < 0.3. 

Hence 



B(t -» a (1450) W)I / T ) <lx 10- 



(25) 



If the ao(1450) is a radial excitation, which is the 
case if the ao(980) is the ud ground state, then 
B(t — > ao(1450) ( - ;rj? ,- ) ^ T ) should be suppressed by an ad- 
ditional wave-function overlap factor. 

Next, we look at the contribution of the eto(980) to 
r — > wr]'i> T , which can be extracted from the relation 



B(r -> ^o(980) (m/) ) = Vfr -> wp(980y) 
B(t -> 1/00(980)^)) ~ U(r - i/ O0 (980) (t)j) ) 

where 



R v = 



7WM980) mf) 



A4(a (980) -> ttt?) 



(26) 



(27) 



is the square of the ratio between the relevant hadronic- 
decay matrix elements. We assume that equals the 
corresponding ratio of ao(1450)-decay matrix elements, 
and is hence obtained from 



2?(a (1450) - mf) Pv 



£(a (1450) -» tttj) p„ 



^, (28) 



where px is the ao(1450)-rest- frame momentum of the 
products of the decay ao(1450) — > 7rX. Given the ~ 50% 
error [l3l | on the ratio of branching fractions appearing 
in Eq. (f28f and the uncertainty on the ao(1450) width, 
R 7 ^ comes out in the range [0.25, 1.25]. The ratio of the 
phase-space integrals in Eq. (|26|) is 0.06, with some de- 
pendence on what one takes for the eto(980) width. Using 



the range for B from Eq. ([2]), we obtain 

B(t -> a (980) (7n/) Jv) w [0.2 to 1.2] x 10" 6 . (29) 



IV. CONCLUSIONS 



Combining Eqs. H2), (J2SJ, and we obtain 

the branching fraction limit 



B(r -> 7rr/V T ) < 1.4 x 10" 



(30) 



in no conflict with the experimental upper limit, Eq. (|3j). 
which is about five times greater. Our result is dominated 
by the ao(980) contribution, assuming it is a ud state. 

The experimental limit was obtained with only a third 
of the currently available BABAR and Belle data sets, and 
with the rj reconstructed only in the 77 final state. There- 
fore, an improvement in the limit can be expected from 
the current generation of B factories, but probably not 
to the level of Eq. ([30]) . By contrast, a Super B fac- 
tory [l4l ] , with two orders of magnitude more luminosity, 
will be able to use B and B' to investigate the nature of 
the ao (980) and to search for new interactions mediated 
by heavy scalars [l|. 
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